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The  surface  structure  of  body-centered-cubic  Fe|210j  has  been  determined  by  low-energy- 
electron-diffraction  intensity  analysis.  Substantial  relaxations  from  bulklike  structure  in  directions 
both  perpendicular  and  parallel  to  the  surface  were  found.  The  |210|  surface  is  the  most  open  of 
the  six  Fe  surfaces  that  have  been  studied,  shows  the  largest  perpendicular  relaxation,  involves  the 
most  structural  parameters  (8),  and  has  been  determined  to  the  greatest  depth  (down  to  the  fifth 
layer).  The  observed  structure  is  as  follows;  spacing  between  first  and  second  layer 
d,2=0.50±0.03  A  (bulk  spacing  is  0.641  A),  in  similar  notation  d2j=0.57±0.03  A, 

1^34=0. 75±0.03  A,  4/45  =0.61  ±0.03  A;  change  of  r^stry  shift  between  first  and  second  layer 
A<Ii2=0.  14±0.05  A,  in  similar  notation  An23=0.03±0.05  A,  Aa34=0.00±0.05  A,  Aa4s 
=0.08±0.05  A. 


I.  INTRODUCTION 

Much  recent  work  by  lov^mergy-electron  diffraction 
(LEED)  intensity  analysis  has  established  accurately  and 
reliably  the  relaxations  of  the  surface  layers  of  a  number 
of  clean  metal  surfacesT^^  where  the  term  relaxation 
refers  to  rigid  translations  of  the  surface  layers  from  their 
bulk  positions  without  change  of  the  unit  cell  of  the  sur¬ 
face  mesh.  Our  work^on  higheij^ndex,  less  symmetri¬ 
cal,  more  open  surfaces  of  bcc  iron,  i.e.,  Fej211), 
Fe{310),  and  Fe(210!,  has  demonstrated  the  occurrence 
of  parallel  as  well  as  perpendicular  relaxations,  both  of 
which  have  damped  oscillatory  magnitudes,  extend  deep 


here  provide  valuable  information  for  further  develop¬ 
ment  of  the  fundamental  theory  of  metal  surfaces.  In  ad¬ 
dition,  we  note  that  structural  changes  of  the  magnitude 
found  here  can  no  longer  be  regarded  as  small  and  it  is 
clearly  necessary  to  fix  the  atom  positions  in  the  surface 
layers  before  quantitative  theories  of  surface  properties, 
e.g.,  vibrational  models  and  electronic  states,  are  possible. 

The  Fe|210|  surface  was  found  to  have  relaxation  ef¬ 
fects  down  to  at  least  the  fifth  layer  including  displace¬ 
ments  both  parallel  and  perpendicular  to  the  surface.  Sec¬ 
tion  II  gives  experimental  details;  Secs.  Ill  and  IV 
describe  the  calculations  and  structure  analysis,  respec¬ 
tively. 


into  the  surface,  and  increase  in  magnituc^e  with  surface  .  ^ 
openness  or  roughness.  ^ ^ 


The  present  work  gives  details  of  the  analysis  of  the 
Fe{210j  surface,  which  is  the  most  open  of  the  six  sur¬ 
faces  of  Fe  that  have  been  studied,  involves  the  largest 
number  of  structural  parameters,  and  is  found  to  have  the 
largest  perpendicular  relaxation  (22%  contraction)  of  any 
of  these  surfaces.  This  large  perpendicular  relaxation  is  in 
qualitative  agreement  with  the  idea  of  a  relaxation  driven 
by  a  surface  Madelung  foice  arising  from  the  flattening  of 
the  electron  density  at  the  surface,  hence  a  Madelung 
force  which  is  larger  for  more  open  surfaces.*  The 
Madelung  force  model  has  been  developed’  by  making 
simple  assumptions  about  the  electron  density  and  the 
electron  response  and  then  allowing  the  lattice  to  relax  so 
as  to  minimize  the  total  energy.  These  models  have  led  to 
interesting  qualitative  features  of  the  relaxation  to  com¬ 
pare  with  measured  relaxations,  as  will  be  noted  later. 
However  deduction  of  the  ground  state  from  first  princi¬ 
ples  by  relaxation  to  self-consistency  of  both  the  electron 
density  and  the  nuclear  positions  has  not  yet  been 
achieved,  even  for  a  simple  metal;  a  transition  metal  like 
Fe  would  be  still  more  difficult.  Hence  the  results  found 


II.  EXPERIMENTAL 

The  methods  for  preparing  an  oriented  single-crystal 
Fel210l  sample  and  the  procedure  for  cleaning  the  crystal 
in  vacuo  were  the  same  as  for  Fe{310).’  Thirty-five 
LEED  intensity- versus-energy  spectra  (31  nondegenerate) 
were  collected  with  a  spot  photometer:  sixteen  spectra  at 
normal  incidence  (01,  1(^  11,  11,  10,  20,  11,11,  20,  30, 
2I,  21,  02,  30,  21,  a^nd  2J),  ten  spectra  at  0=7. 4°, (^=90° 
(00,  10,  10,  11,  11,  1  1,  2  l,_21,_20,_and  30^Md  nine  spec¬ 
tra  at  0=  13.2",<^=90“  (00,  10,  20,  30,  31,12,  21,  ll,  and 
32).  All  data  were  collected  with  the  crystal  at  room  tem- 
ture  and  the  angles  and  beam  indices  follow  the  conven¬ 
tion  of  Zanazzi  et  al.*  A  schematic  drawing  of  the 
LEED  pattern  is  shown  in  Fig.  1 . 

III.  CALCULATIONS 

LEED  intensity  spectra  were  calculated  with  the  com¬ 
puter  program  change.’  *®  The ^ very  small  interlayer 
spacing  of  Fej  210)  ( dbun  =0.641  A)  and  the  correspond¬ 
ingly  large  surface  unit  cell  (2.87  A X  6.41  A)  necessitated 
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FIG.  1.  Schematic  LEED  pattern  of  FefllO)  for  0=0*, 
^=90*. 

the  use  of  up  to  123  plane  waves,  or  beams,  to  represent 
the  wave  function  accurately  for  calculating  interlayer 
scattering.  Up  to  five  planes  of  atoms  were  bunched  to¬ 
gether  into  composite  layers,  among  which  the  scattering 
was  calculated  in  the  spherical  wave  basis,  rather  than  the 
beam  basis,  and  hence  was  not  restricted  to  widely  spaced 
layers.  The  Fe  potential  was  the  same  as  used  for  previ¬ 
ous  iron  work."  Eight  phase  shifts  were  used  to  represent 
ion  core  scattering  and  the  m^n  vibrational  amplitude 
was  taken  as  (<m^))'''^=0.115  A.  The  complex  inner  po¬ 
tential  Vi  —  Vq—IP  was  taken  as  energy  independent,  P 
was  set  up  to  4  eV,  and  Vq  was  left  to  be  determined  from 
the  intensity  analysis  with  the  initial  value  fixed  at  — 11.5 
eV. 


TABLE  I.  Surface  lattice  vector  aj,  and  82,  parallel  com¬ 
ponent  Cp,r  of  the  interlayer  translation  vector,  interlayer  dis¬ 
tance  d,  and  mirror  plane  m  for  the  bulklike  Fef210j  surface. 
1  and  ^  are  directed  along  <210)  and  <  100),  respectively,  as  la¬ 
beled  in  the  top  view  of  Fig.  2;  t  is  directed  along  the  inward 
surface  normal,  a  =2.866  A. 
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IV.  STRUCTURE  ANALYSIS 

The  bulklike  Fe|  210j  surface  is  stepped  and  consists  of 
1 1  lOj  terraces  four  atomic  rows  wide  joined  by  monatom¬ 
ic  |110j  steps.  Top-  and  side-view  diagrams  of  the  sur¬ 
face  are  shown  in  Fig.  2;  Table  I  gives  the  structural  pa¬ 
rameters.  The  unit  cell  is  a  rectangle  with  lattice  vretors 
<1  directed  along  a  (210)  axis  and  of  length  6.41  A  and 
02  along  a  <  100)  axis  and  of  length  2.87  A.  The  inter¬ 
layer  translation  vector  c  is  in  a  (closed-packed)  (111) 
direction  with  |c|  =2.48  A.  The  interlayer  spacing  is 
0.641  A  and  the  bulklike  registry  of  successive  layers  is 
described  by  the  relation;  Cp„=3ai/10-|-a2/2,  where  Cp,, 
is  the  projection  of  c  onto  the  j210)  plane.  The  surface  is 
veiy  open,  with  a  packing  fraction  of  0.2634  (the  most 


FIG.  2.  Top  and  side  views  of  the  undistorted  Fe|210j  surface.  Circles  of  equal  thickn«s  represent  coplanar  atoms,  for  the  top 
view  circles  of  decreasing  thickness  indicate  progressively  deeper  layers.  All  distances  are  in  A. 
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FIG.  3.  Schematic  of  the  top  four  layers  of  the  Fe|210{  sur¬ 
face.  Circles,  first  layer;  triangles,  second  layer;  squares,  third 
layer;  hexagons,  fourth  layer. 


densely  packed  bcc  surface,  f  llOj,  has  a  packing  fraction 
of  0.8330).  Six  layers  are  visible  at  the  surface  in  a  touch¬ 
ing  hard-sphere  model  and  fourth-layer  atoms  are  at  a 
nearest-neighbor  distance  (2.48  A)  from  their  closest 
neighbors  in  the  top  layer  (see  Fig.  3).  Hence,  there  are 
many  nonequivalent  atoms  which  can  be  legitimately 
described  as  surface  atoms,  a  fact  which  makes  the  (210| 
face  both  interesting  and  complex.  For  example, 
numerous  different  adsorption  sites  exist  on  the  surface 
for  chemisorbed  monolayers.  We  have  performed  experi¬ 
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FIG.  4.  Variation  of  the  r  factor  for  the  0=0°  data  set  as  a 
function  of  dn  and  du  with  dj4=0.748  A  (0.735A),  d45 
=0.606  A  (0.621  A),  0,2  =  2.056  A,  023  =  1.958  A  (1.921  A), 
034=1.942  A  (1.979  A),  045=1.995  A,  and  Ko=-11.5eV. 


ments  studying  oxygen  and  sulfur  structures.  Both  oxy¬ 
gen  and  sulfur  form  (2x1)  overlayers  but  as  yet  the  struc¬ 
tures  of  these  systems  remain  unsolved. 

The  structural  models  considered  included  variations  of 
the  first  four  interlayer  spacings,  d|2,  <^23,  ^34,  and 
1/45  (dbuiy  =0.641  A),  and  the  first  four  registry  shifts, 
0(2,  a2u  034,  and  043  (the  registry  shift  between  succes¬ 
sive  layers  /  and  J  is  defined  by  Cp,,  =a2/2-t-ay'i,  the 
bulk  value  of  Ofj  is  1.923  A). 

The  first  series  of  calculations  involved  independent 
variations  of  the  first  three  interlayer  spacings  and  the 
top-layer  registry,  by  relatively  large  amounts,  in  order  to 


TABLE  II.  Results  of  calculations  for  6=0’  of  various  structural  models:  d,:,  d^,  and  d-^  are  the 
first  three  interlayer  spacings,  0,2  and  au  are  the  first  two  interlayer  registry  parameters  (see  text).  The 
agreement  between  theory  experiment  is  given  by  the  reliability  factor  r.  All  distances  are  in  A. 


Run 

d,2 

d23 

d34 

O12 

023 

034 

r 

Kbulk) 

0.641 

0.641 

0.641 

1.923 

1.923 

1.923 

0.293 

2 

0.561 

0.641 

0.641 

2.188 

1.923 

1.923 

0.216 

3 

0.587 

0.614 

0.641 

2.320 

1.923 

1.923 

0.236 

4 

0.561 

0.614 

0.641 

2.452 

1.923 

1.923 

0.249 

5 

0.561 

0.667 

0.641 

2.188 

1.923 

1.923 

0.226 

6 

0.561 

0.640 

0.667 

2.320 

1.923 

1.923 

0.235 

7 

0.534 

0.587 

0.641 

2.188 

1.923 

1.923 

0.208 

8 

0.534 

0.587 

0.641 

2.055 

1.923 

1.923 

0.185 

9 

0.561 

0.587 

0.641 

2.053 

1.923 

1.923 

0.194 

10 

0.561 

0.693 

0.641 

2.188 

1.923 

1.923 

0.270 

11 

0.534 

0.640 

0.641 

2.188 

1.923 

1.923 

0.187 

12 

0.534 

0.693 

0.641 

2.188 

1.923 

1.923 

0.270 

13 

0.540 

0.667 

0.640 

2.188 

1.923 

1.923 

0.179 

14 

0.518 

0.667 

0.640 

2.188 

1.923 

1.923 

0.164 

15 

0.518 

0.614 

0.640 

2.188 

1.923 

1.923 

0.154 

16 

0.518 

0.614 

0.667 

2.188 

1.923 

1.923 

0.144 

17 

0.518 

0.598 

0.693 

2.188 

1.923 

1.923 

0.142 

18 

0.508 

0.598 

0.693 

2.10s 

1.923 

1.923 

0.138 

19 

0.4S7 

0.598 

0.693 

2.108 

1.923 

1.923 

0.132 

20 

0.503 

0.598 

0.730 

2.108 

1.923 

1.923 

0.127 

21 

0.503 

0.598 

0.730 

2.108 

1.976 

1.923 

0.127 

22 

0.503 

0.598 

0.730 

2.108 

2.029 

1.923 

0.128 

23 

0.503 

0.598 

0.704 

2.055 

1.997 

1.923 

0.124 

24 

0.496 

0.581 

0.708 

2.056 

1.995 

1.923 

0.118 

25 

0.496 

0.581 

0.708 

2.056 

1.995 

1.976 

0.115 

1932 


J.  SOKOLOV,  F.  JONA,  AND  P.  M.  MARCUS 


31 


TABLE  III.  r  factors  for  the  bulklike  and  fully  relaxed 
models  for  data  collected  at  three  angles  of  incidence. 


Data  set 

r 

Bulklike 

factors 

Relaxed 

Energy  range  (eV) 

6=or 

0.293 

0.103 

1958 

fl-7  4*1 

4=90*  1 

0.207 

0.113 

928 

fl- 13.2*1 

4=90*  1 

0.347 

0.125 

716 

Total 

0.282 

0.110 

3602 

get  a  general  idea  of  what  the  optimum  structure  might 
be.  Agreement  between  theory  and  experimental  spectra 
was  measured  using  the  numerical  reliability  factor  of 
Zanazzi  and  Jona.'^  For  each  calculation,  the  value  of  the 
nonstructural  parameter  Vq,  the  real  part  of  the  itmer  po¬ 
tential,  was  allowed  to  vary  independently  from  —7.3  to 
— 15.5  eV;  the  best  value  was  always  close  to  — 11.5  eV. 
A  sampling  of  the  first  set  of  calculations  (about  y  of  the 
total  number  of  calculations)  is  shown  in  Table  II  labeled 
with  run  numbers  1—20.  Next,  preliminary  calculations 
were  made  varying  the  second  and  third  registry  shifts, 
023  and  <*34  (runs  21—25  of  Table  II  are  a  sampling).  At 
this  stage,  the  calculations  indicated  that  the  b«t  models 
were  those  with  contracted  (2-20%),  ^23  contract«l 
(^10%),  d}4  increased  (2£ll%),  ai2,  023,  and  a-^  in¬ 
creased  by  approximately  7%,  4%,  and  3%,  respectively, 
compared  to  the  bulk  values.  The  next  series  of  calcula¬ 
tions  fixed  all  but  two  of  the  structural  parameters  (set 
near  to  the  corresponding  values  found  above)  and  then 
allowed  the  remaining  two  parameters  to  vary  indepen¬ 
dently  until  a  two-dimensional  minimum  of  the  reliability 
factor  was  located.  The  best  values  of  the  two  variables 
were  determined  by  fitting  the  r  (reliability)  factors  to  a 
quadratic  function  (elliptic  paraboloid)  near  the 
minimum.  With  the  best  values  thus  obtained,  a  new  pair 
of  variables  was  chosen  and  the  process  repeated.  Over¬ 
lapping  of  pairs  was  considered,  for  example 
(d}4,d4s),  and  (dijtduY,  when  the  consistency  among  op¬ 
timum  values  determined  as  members  of  different  pairs 
was  as  good  or  better  than  the  estimated  experimental  un¬ 
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FIG.  5.  Schematic  of  relaxed  Fe|210|  surface  with  optim¬ 
ized  structural  parameters.  All  distances  are  in  A. 


certainties  (for  a  discussion  of  uncertainties  in  LEED,  see 
Sokolov  et  al.*),  the  analysis  was  considered  finished.  An 
example  of  the  sensitivity  of  the  r  factor  to  the  parameter 
values  is  given  in  Fig.  4  for  the  pair  (di2,d2i).  For  the 
variables  whose  best  values  were  determined  as  members 
of  different  pairs  (023.  <*34.  <^23>  <^34)  average  of 

the  two  “best”  values  was  taken  for  the  final  results  quot¬ 
ed  below.  In  all  ca^,  the  difference  between  the  two  best 
values  were  0.01  A  or  less.  The  results  of  the  analysis 
(with  estimated  errors  following  Sokolov^  et  al.*)  were  as 
follows:  (bulk  interlayer  spacing  0.641  A,  bulk  layer-to- 
layer  registry  shift  1.923  A) 

di2=0.50±0.03  A  (22.0±4.7%  contraction), 
d23=0.57±0.03  A  (11.1  ±4.7%  contraction), 
d34=0.75±0.03  A  (17.0±4.7%  expansion), 

^45  =0.61  ±0.03  A  (4.8±4.7%  expansion), 
ai2  =  2.06±0.05  A  (7.1  ±1.6%  increase), 

023=  1.95±0.05  A  (1.4±2.6%  increase), 
o„=1.92±0.05  A  (0±2.6%) , 
a45=2.00±0.05  A  (4.0±2.6%  increase), 
Fo=-11.5±1.0eV  , 

'■min=0103  (r  factor  of  Zanazzi  and  Jona'^or  16 
beams  at  6=Cf). 

The  percentage  changes  are  with  respect  to  the  bulk 
values.  Figure  5  shows  a  schematic  representation  of  the 
optimized  structure. 

Calculations  were  then  made  with  the  above  optimum 
values  for  the  data  sets  taken  at  off-normal  incidence.  Be¬ 
cause  of  the  extreme  length  of  the  calculations  (about  4  h 
of  CPU  time  on  an  IBM  3081  computer),  it  was  not  possi¬ 
ble  to  repeat  the  analysis  for  the  other  data  sets. 

Table  III  gives  the  r  factors  for  both  the  bulklike  and  ful¬ 
ly  relaxed  models  for  all  35  /-vs-  V  spectra  and  Figs.  6—8 
show  the  plots  of  the  corresponding  curves. 

It  should  be  emphasized  at  this  point  that  this  structur¬ 
al  problem,  involving  eight  structural  parameters,  is  not¬ 
ably  more  difficult  than  any  other  attempted  so  far  for 
metal  surfaces.  Thus,  the  optimization  process  carried 
out  in  this  work  is  not  exhaustive,  and  further  refinement 
would  be  expected  to  improve  the  agreement  betw^ 
theory  and  experiment ^parti<'nlarly  for  spectra  like  the  2  1 
beam  at  6=(f  and  the  30  beam  at  6=7.4°,  ^=90*). 

To  summarize,  the  loosely  packed  Fe{210|  surface  was 
found  to  have  large  relaxations,  both  normal  atid  parallel 
to  the  surface  plane.  The  pattern  of  contraction,  contrac¬ 
tion,  expansion,  and  contraction  for  di2,  d23,  d^,  and 
d4j,  respectively,  is  in  agreement  with  the  form  predicted 
for  bcc  Naj210j.’  However,  our  results  for  the  parallel 
relaxations  (all  parallel  motions  in  the  same  direction)  do 
not  agree  with  the  results  of  Barnett  et  al.''  The  calcula¬ 
tions  did  not  include  the  effects  of  screening  of  the  ions 
by  valence  electrons  and  so  it  is  an  opoi  question  as  to 
whether  a  more  refined  theory  would  give  results  closer  to 
the  observed  Fe{210j  structure  or  differences  between  the 
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FIG.  7.  Same  as  for  Fig.  6  but  for  0=7.4°,  ^  =  90°. 


50  so  70  90  110  150  ISO  170  190  50  SO  70  90  110  ISO  ISO  170  1*0  50  SO  70  90  110  ISO  ISO  170  190 

ENEROY  UV)  ENERGY  <wl  ENERGY  IW) 

FIG.  8.  Same  as  for  Fig.  6  but  for  0=  13.2°,  4=90°. 


electronic  structure  of  Na  (a  simple  metal)  and  Fe  (a  tran¬ 
sition  metal)  are  important. 
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